Bayesian analysis of Friedmannless cosmologies 
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Assuming only a homogeneous and isotropic universe and using both the 'Gold' Supernova Type 
la sample of Riess et al. and the results from the Supernova Legacy Survey, we calculate the Bayesian 
evidence of a range of different parameterizations of the deceleration parameter. We consider both 
spatially flat and curved models. Our results show that although there is strong evidence in the 
data for an accelerating universe, there is little evidence that the deceleration parameter varies with 
redshift. 



I. INTRODUCTION 

The accelerating expansion of the universe is probably 
the most important discovery in cosmology in the last 
decade. Direct evidence for cosmic acceleration comes 
from the Hubble diagram with supernovae of type la 
(SNIa) as standard (or standardizable) candles 0, 0, Q , 
and indirect support comes from e.g. observations of the 
large-scale distribution of matter [J, [j| combined with 
measurements of temperature anisotropics in the cosmic 
microwave background (CMB) radiation [6|. 

While the most straightforward interpretation of the 
acceleration is that the energy density of the universe is 
presently dominated by the cosmological constant, the 
problem of understanding why it is so small compared 
with the natural energy scale of quantum gravity has led 
to proposals of a jungle of alternatives. The proposals can 
be roughly divided into two classes: modifications of the 
right-hand side or the left-hand side of the Einstein equa- 
tions. The first case corresponds to changing the theory 
of gravity from standard general relativity (GR), while 
in the second case some kind of negative-pressure fluid is 
added to the energy-momentum tensor. At present none 
of these possibilities can be ruled out. The literature is 
vast, but useful reviews are found in references 

In this paper we ask the question of what one can learn 
from the SNIa data about the kinematics of the universe. 
This sidesteps the issue of the cause of the acceleration 
and leads to fairly model-independent conclusions since 
we only assume that the universe is homogeneous and 
isotropic and is described by a metric. We think a study 
like this is timely since sometimes one sees in the lit- 
erature wide-ranging conclusions drawn on the basis of 
fitting specific models to SNIa data. For example, the 
recent discussion of 'phantom energy' ^(J is based on 
the fact that fits with a dark energy component with a 
constant equation of state parameter w prefer values less 
than —1. We think it is important to bear in mind what 
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exactly the SNIa data is telling us about the expansion 
of the universe when assessing the merit of claims like 
this. 

Our study is similar in spirit to that of a recent pa- 
per by Shapiro and Turner |7(, but our analysis is based 
on the concept of the Bayesian evidence. Also, we an- 
alyze both the 'Gold' SNIa sample of Riess et al. and 
the more recent sample from the Supernova Legacy Sur- 
vey (SNLS), and show that they lead to slightly different 
conclusions about the expansion history of the universe, 
although none of the data sets can provide strong con- 
straints. 



A. Friedmannless Cosmology 

Assuming homogeneity and isotropy on large scales, 
the requirement for a maximally symmetric subspace 
leads us to the Friedmann-Robertson- Walker (FRW) 
metric: 



ds 2 



dt 2 -a 2 (t) 



dr 



1 - kr 2 



r 2 dfl 2 



(1) 



where a(t) is the scale factor. Note that we have not 
assumed spatial flatness, k = 0, in order to keep the 
discussion as general as possible. 

The deceleration parameter is defined as 
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where we have written q in terms of the Hubble param- 
eter and ' = d/dz. Similarly, the jerk is given by: 
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Naturally, one may consider expansions of q(z), j(z) or 
even higher order parameters. However, the order of the 
parameter, or the number of derivatives, is related to 
the number of free parameters in H(z). This is clearly 
demonstrated by considering for example models with 
constant deceleration parameter, q = qo, or jerk, j = jo. 
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In this case equations |J2J and are straightforwardly 
solvable with solutions of the form 



H 2 (z) 
H 2 (z) 
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where 



*i,2 = l± \/2(-io-l) + 7- 



(4) 
(5) 
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Hence, expanding the jerk can be viewed as requiring 
more parameters, or cosmic fluids, than expanding the 
deceleration parameter. Guided by these considerations 
we consider expansions of q(z). This furthermore allows 
us to compare our results with those of who reach 
their conclusions by a different method. 

In a standard matter only universe, Einstein-de Sit- 
ter (EdS) model, the deceleration parameter is simply 
q = 1/2 and jerk j = — 1. A ACDM universe with 
H 2 = i?o(fi m (l + z) 3 + Q\) on the other hand has a 
non-constant deceleration parameter, 



q(z) = - 



n A /n m -^(i + z ) 3 
n A /n m + (i + z) 3 



(7) 



and again a constant jerk j = — 1. This also demonstrates 
how mapping from cosmological models to cosmological 
parameters is not unique, eg. both EdS and ACDM mod- 
els have the same jerk further supporting the use of q as 
an expansion parameter. 

We make minimal or no assumptions about the matter 
content of the universe. Considering the deceleration pa- 
rameter we can avoid the question of the gravitational 
theory, whether standard Einstein's equations hold or 
not, ie. we make no assumptions on the connection be- 
tween the energy density of matter and the evolution of 
the universe. Similar constructions have been considered 
previously: Dvali and Turner ^jj considered adding a 
correction to the Friedmann equation, Freese & Lewis 
|l2j considered generalized Friedmann equations in their 
Cardassian models and generalizations in a similar spirit 
were also considered in a previous paper 

111 

In addition to considering straightforward linear ex- 
pansions of the deceleration parameter, q(z) = ^qiZ 1 , 
we also study more physically guided expansions, q(z) = 
9i(l+z)~ 31 . Such expansions are motivated by the fact 
that if we assume that the energy density is conserved in 
the expanding universe, cold dark matter density is also 
conserved and hence p m ~ a~ 3 ~ (1 + z) 3 . One can 
view this requirement simply as arising from conserva- 
tion of the number of the particles and not as a property 
of the Einstein's equations. Energy conservation for ex- 
ample holds in a large class of modified gravity models, so 
called f(R) models both within the metric and Palatini 
approach [jjj- As example, consider the ACDM model, 
which when expanded in powers of (1 + z)~ 3 gives 



i_3^V (1 + 0) -3 
2 2 m 



The models we consider along with the very conserva- 
tive priors are shown in table Q] 



Model 
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Mi 

M 2 
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M 4 

M, 



Non-flat 



M 6 



M 7 



q = qo 

q = qo + qiz 

q = q Q + q x z + q 2 z 2 

q = qo + qi/(l + z) 3 



q = <?o 



q = 



+ qi /(l + z) 3 + q 2 /(l + zf 

qo, z < zt 
qi, z > z t 



q = qo 



q = qo + qiz 



Priors 



TABLE I: Model parameters and priors 



B. Bayesian Evidence 

The Bayesian evidence (BE), or E(AI), is defined as 
the probability of the data D given the model M with a 
set of parameters 9, 

E{M) = P{D\M) = J d8P(D\6,M)P(0\M), (9) 

where P(0\M) is the prior on the set of parameters, nor- 
malized to unity. We follow here the common assumption 
that the parameter priors are top hat and hence © can 
be written as 



E(M) 



Jd6P{D\6,M) 



(10) 



In this paper we approximate the probability by 
P(D\9) « cxp(— x 2 (0)/2). Since we are interested in 
comparing models, we leave the normalization arbitrary. 

In addition to BE, we will at times refer to an unnor- 
malized Bayesian evidence, which is simply 



E(M) = / d6P(D\6,M). 



(11) 



3 



In order to compare models, we utilize the Bayes Fac- 
tor which is defined as the ratio of evidences of two mod- 
els and M,-: 



where simifc(x) is the commonly used function 



Bi 



E{Mj 
E(M j 



(12) 



If Bij > 1, model Mj is preferred over Mj , given the data. 
In similar spirit as in the frequentist approach, where 
one compares the x 2 s of given models and concludes how 
significantly a model fits the data better, in the Bayesian 
Evidence framework one can assign significance to the 
difference between models by using the Jeffreys Scale Q 
as shown in Table CU Note that Ba is assumed to be 



Log of Bayes Factor 


Significance 


In Bij < 1 


Not significant 


1 < \nB t j < 2.5 


Substantial 


2.5 < lnB^ < 5 


Strong 


5 < In 


Decisive 



TABLE II: Jeffreys Scale 



greater than one ie. model Mj has larger evidence than 
model Mj. 

An approximation to the Bayesian Evidence is the 
Bayesian Information Criterion (BIC) [l5lll^ |. BIC also 
penalizes for extra parameters by reducing the likelihood, 



BIC = -2£ + k InN, 



(13) 



where k is the number of parameters and N is the number 
of data points used in the fit. Approximating the likeli- 
hood by x 2 ; it is evident how in the frequentist language 
the fit to a model with more parameters needs to be sub- 
stantially better in order to justify the extended model. 
For example, for a SNIa dataset with 157 supernovae, a 
two parameter model must have a significantly better fit 
to the data, A\ 2 ~ 10, than a single parameter model in 
order to justify the more extended model (ABIC w 5). 



II. METHODS 

In this work we use data from two SNIa compilations, 
the Riess Gold set Q with 157 supernovae up to z = 1.76 
and the more recent SNIa data released by the Super- 
nova Legacy Survey (SNLS) i containing 115 super- 
novae with z < 1. 

Within the Friedmannless Cosmology approach, the 
geometry of the universe becomes relevant at this stage 
since the luminosity distance depends on the spatial cur- 
vature k. The dimcnsionlcss luminosity distance is 



d L (z) = (1 + z)sinn fc 



du 
H(uj 



(14) 



simifc (x) 



sin(x), k = +1 
x, k = 
sinh(x), k = — 1 



(15) 



The two surveys use somewhat different methods for 
calculating the \ 2 of a model (see eg. 0, El El f° r a 
more detailed description). For the Gold set, the x 2 of a 
model M with a set of parameters {on} is 

2/r i \ (AW» - 51o g 10 ci L (z t ; {on}) - Ho) 2 

x (\ a ih Mo j - 
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(16) 

where is the Hubble free luminosity distance. The 
Hubble parameter and corrections to the absolute mag- 
nitude arc hidden in /iq- 

The analysis for the SNLS dataset is somewhat more 
complicated due to the appearance of two extra nuisance 
parameters a and /3: 



E 



X 2 ({ai},M + no,a,/3) 
(Vobs,i - 5 log 10 d L (zi; {ai}) - ^ ) 2 



(17) 



where now er^ = cri(a,j3). The computation of \ 2 us- 
ing the SNLS is hence computationally more intensive as 
there are three nuisance parameters to marginalized over 
(again, we refer to 0] for a detailed description of the 
procedure). 

An often utilized approximation in marginalizing over 
nuisance parameters is to marginalize by maximizing the 
likelihood. For example, instead of marginalizing over /j,q 
by integrating over fiQ with a given prior, a commonly 
used approximation when using the Gold set in study- 
ing cosmological models is to maximize the likelihood by 
minimizing \ 2 with respect to ^o- This is well motivated 
as can be seen by first expanding l|16|) . which gives 



X 2 ({ai},Ho) = ci - 2c 2 ^q + c 3 ^q, 
where Ci are independent of /j,q : 

E (Hobs,i - 51og 10 d L (zi- {aj})) 7 



(18) 



C2 = E 

i 
n 

<3 =£ 



fJ-obs,i - 5 log 10 dh(zi\ {ai}) 



(19) 



The minimization over /xo is trivial, fiQ — C2/C3. Hence, 
instead of finding the minimum of X 2 ({ a i}i Mo) w e can 
minimize X 2 ({ a i}) = X 2 {{ a i\ i °2 / c z) ■ In calculating the 
confidence contours, however, this method is only ap- 
proximate unless we assume no prior on fio ■ I n the special 
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case with no prior on one can carry out the integral 
explicitly resulting in a constant term that will not have 
an effect in calculating the confidence contours. 

As a check of our method, we have confirmed that the 
confidence contours produced using the Gold sample for 
the two epoch model, M 5 , agree with those presented in 
■ For the SNLS sample, we checked that we reproduced 
the results in Q for flat and non-flat ACDM models. 

In calculating the Bayesian evidence E{M) of a given 
model, integration over the parameters is required. In 
order to assess the validity of the approximation of 
marginalizing over nuisance parameters by maximizing 
the likelihood we calculate E(M) for the Gold set both 
by using explicit integration over /io and by the approx- 
imate method. The parameter priors are chosen to be 
qo £ [-1,1], qi e [-2,2], fi & [42.38,43.89] and the 
comparison is carried out for the constant q and linear q 
model. When marginalization over /j,q is done explicitly 
when calculating BE, we expect the resulting evidence 
to be less than when /i is marginalized by maximizing 
the likelihood at each point. This is due to normalization 
and to the fact that the likelihood at a point in parame- 
ter space is typically less when not explicitly maximized. 
This proves to be the case, eg. for the constant model 
E(Mq) is about 30 times larger when the marginalization 
is carried out approximatively. However, in using BE we 
are interested in comparing models and hence a more 
interesting issue is to compare the two different models 
with both marginalization methods. We find that the 
difference between marginalization schemes is not signif- 
icant when comparing models, £?oi ~ 2.5 in both cases. 
Hence, marginalizing over /iq by maximizing the likeli- 
hood appears to be a reasonable approximation in calcu- 
lating the Bayesian Evidence at least for the parameter 
priors used in this paper. 



Model 


In B(h 


\ 2 Best fit parameters 


M 





182.8 go = 


-0.29 


Ivl 1 


9 n 
z.u 


1(0.4 < 


' qo 


= -0.70 








v 9i 


= 1.5 








' 9o 


= -1.1 


M 2 


1.7 


173.8 { 


<?i 


= 4.5 








>. 92 


= -3.3 






174.0 < 


qo 


= 0.70 










= -1.8 








' 9o 


= -0.30 


M 4 


2.6 


172.7 « 


* 


= 3.2 








, 92 


= -4.9 








' 9o 


= 1.6 


M 5 


1.8 


172.4 < 




= 0.20 










= 0.11 



Non-fiat 



M 6 , 


(* = 


-1) 


3.5 


176.2 


9o 


-0.64 


M 6 , 


(* = 


1) 


< 


191.1 


7 = 


-0.04 


M 7 , 


(fc = 


-1) 


1.9 


176.2 


\ 


9o 


= -0.70 












I 9i 


= 0.20 


M 7 , 


(fc = 


1) 


1.5 


176.0 


\ 


9o 


= -0.70 










I 9i 


= 2.3 



TABLE III: Bayes factor relative to the constant q model, 
minimum \ 2 and best fit parameters for different models 

ima. This further advocates the use of Bayesian Evidence 
as an efficient model selection tool. 



B. The SNLS sample 



III. RESULTS AND DISCUSSION 
A. Gold sample 

The Bayes factor relative to the flat constant q model, 
ln£?oi, minimum % 2 and best fit parameters for the dif- 
ferent models are given in table IIIII Keeping in mind 
the Jeffrey's scale, we see from the table that most mod- 
els are not significantly better than the baseline model 
Mo, which is the constant q flat universe model. Only 
flat models for which we find strong evidence are those 
where we have assumed a more physically motivated ex- 
pansion (note that A In .By = 5 corresponds to odds of 
1/13). Interestingly, the model with the largest evidence 
is a constant q closed universe, which is strongly preferred 
over its flat universe counterpart. 

In calculating the Bayesian Evidence we use a simple 
grid. We have checked that the results are insensitive to 
further refinement of the grid. Note that in contrast the 
best fit parameters maybe sensitive to grid spacing, for 
example due to the existence of nearly degenerate min- 



Carrying out the same calculations for the SNLS sam- 
ple of SNIa gives the results in table IIVI The nuisance 
parameters a, (3, and M + /io were in each case fairly 
tightly constrained by the data, so we carried out the 
marginalization over them by fixing them to their best- 
fit values. The model parameters were integrated over 
with the priors given in table I. An unsatisfactory aspect 
of this approach is that, e.g. the likelihood function for 
the Mi model does not reduce to that of M\ for qi = 0, 
or to that of Mq for qi = q2 = 0. We can get a sense 
of how this affects the results by using the best-fit a, (3 
and M + \i for Mq for all the models. The results for this 
case are given in table Note that the Bayes factors 
do not change much except for models M6(k — — 1) and 
M6(k = +1). 

C. Model comparison 

We can now rank the models based on their Bayes 
factor relative to Mq. Both the Gold sample and the 
SNLS provide substantial evidence against all models 
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Mi 
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' g = -0.32 


M 2 
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gi = -2.01 
k g 2 = 4.64 
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M3 




2.0 
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' go = -0.11 
. gi = -0.47 
' g = -1.06 


M 4 
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111.4 < 


gi = -4.97 
. g 2 = 3.73 
' g = -0.46 


M s 
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111.3 < 


gi = 0.35 
k z t = 0.50 


Non-flat 




M 6 , (k = 


-1) 


3.3 


115.2 


go = -0.24 


M 6 , (k = 


+1) 


3.4 


112.2 


go = -0.65 


M 7 , (k = 


-1) 


1.7 


111.5 ■ 


' go = -0.53 
. gi = 1.53 


M 7 , (k = 
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1.5 


111.6 - 


' go = -0.51 
[ gi = -0.73 



TABLE IV: Bayes factor relative to the constant g model, 
minimum \ 2 an d i>est fit parameters for different models from 
the SNLS sample. 



TABLE V: Bayes factor relative to the constant g model, min- 
imum x 2 and best fit parameters for different models from the 
SNLS sample. The nuisance parameters a, /3 and M + fio were 
marginalized over by fixing them to their best-fit values for 
the M0 model. 



from places 3 up to and including 7, and strong evidence 
against models at ranks 8-10, but they do not agree on 
which models fall into which category. However, the best 
model in both cases has q(z) = constant. It therefore 
seems fair to conclude that there is no significant evi- 
dence in the present supernova data for a transition from 
deceleration to acceleration, and claims to the contrary 
are most likely an artifact of the parameterization used 
in the fit to the data. Comparing the ranking accord- 
ing to the BIC gives very different results for both data 
sets, and thus one should avoid approximations like this 
in Bayesian model selection problems. 

That the two SNIa samples prefer slightly different 
models has already been noted in, e.g., [18j where it was 
shown that the SNLS best-fit model was in concordance 
with the best-fit WMAP model , whereas this was not 
the case for the Gold sample. In our case, the Gold sam- 
ple prefers a closed constant-q model, whereas the SNLS 
sample prefers a flat constant-g model. Due to the more 
ho mog eneous nature of the SNLS sample, the authors 
of [18| prefer this sample. For both cases, however, our 
study shows that there is little model-independent infor- 
mation to be extracted from the data beyond the fact 
that the universe is accelerating. 



Rank 


Gold sample 


SNLS 


Gold sample (BIC) SNLS (BIC) 


1 


M 6 (k = +1) 


Mo 


M 6 (fc= -1) 


Me(fc = +1) 


2 


Mo 


Mi 


M 3 


Mo 


3 


M 7 (k = +1) 


M 4 


Mi 


M 6 (k = -1) 


4 


M 2 


Mr(k = +1) 


M 7 (fe = +1) 


Mr(k = +1) 


5 


M 5 


Mi 


M r (k = -1) 


Mr(k = -1) 


6 


M 7 (k = -1) 


M 7 (k = -1) 


Ms 


Mi 


7 


Mi 


M 3 


M ,M 4 


M 3 


8 


M 4 


Me{k= -1) 


M 2 


M 2 ,M 5 


9 


M 3 


M 6 (fc = +1) 


M 6 (fc = +1) 


M 4 


10 


M 6 (fc = -1) 


M 5 







TABLE VI: The models ranked according to Bo; for the Gold 
sample and the SNLS sample, and ranked according to the 
Bayesian Information Criterion. 



D. Redshift of transition from deceleration to 
acceleration 

Riess et al. used the 'Gold' sample to constrain the 
redshift of transition from deceleration to acceleration, 
finding a transition at z t = 0.46 ± 13 0]. This result 
was obtained using a linear expansion of q(z), the model 
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we have labelled Mi. As pointed out by Shapiro and 
Turner []j , this model always has a transition redshift at 
z t = —qo/qi, provided go and gi have opposite signs, so 
one should not assign too much significance to it. In par- 
ticular, the best-fit values of go and gi for a given data 
set might easily give a value of Zt outside of the redshift 
range of the SNIa sample used, and Zt is then clearly 
unphysical. We find this to be the situation with the 
SNLS sample. The best-fit values of go and gi in table 
El give z t fa 2.0, while the SNLS sample only extends 
out to z ~ 1. That the linear expansion is not a good 
way of constraining zt is also evident from the likelihood 
contours in the go-gi plane, shown in figure ^ The data 
are consistent with a wide range of values of zt , even neg- 
ative values are allowed. A better way of looking for a 




FIG. 1: Likelihood contours in the qo-qi plane for the fit of 
model Mi to the SNLS data. The straight lines are lines of 
constant Zt, with Zt varying between —2 and 2. 

transition from deceleration to acceleration is by way of 
our model M 5 : parametrizing q(z) as a piecewise con- 
stant in two bins with the redshift of transition as an 
additional parameter in the fit. However, since neither 
the SNLS nor the 'Gold' sample favour this model, it is 
at the moment not possible to say anything about when 
(or, indeed, if) the Universe went from deceleration to 
acceleration. 



to normalization. In other words Bayesian Evidence pe- 
nalizes a model for large parameter priors, which sounds 
reasonable since then model with a large uncertainty in 
parameter space will typically be classified as less signif- 
icant as the same model but with tighter priors. 

As an example, consider the Bayesian Evidence for the 
constant g model depicted in Fig. Ufa). The (red) solid 
line is the BE as defined in Eq. 10 and the (green) 
dashed line is the unnormalized Bayesian Evidence, Eq. 
(|lip. The prior volume in this case is just the differ- 
ence between the maximum and minimum go value cen- 
tered around go = 0, ie. the horizontal axis represents 
the size of the prior. For this model, the maximum 
likelihood point, or minimum \ 2 = 182.8 was found at 
go = —0.29. From the figure we see that the Bayesian Ev- 
idence first grows quickly with increasing prior volume, 
V p , then reaches a maximum and after that begins to de- 
crease. The unnormalized BE behaves similarly until the 
maximum value at which saturates to a constant. This 
behavior is simple to understand: when the prior is too 
small so that is does not encompass the point of maxi- 
mum likelihood, the BE is relatively small. The evidence 
then grows with the prior volume since the likelihood 
increases as more likely points are included within the 
prior. As the prior becomes large enough to encapsu- 
late the volume where the likelihood is concentrated, the 
normalization factor becomes more important and the 
evidence begins to decrease. As the unnormalized likeli- 
hood does not penalize for the size of the prior volume, 
it simply grows until all of the significant likelihood is 
within the integral and then asymptotes to a constant. 



100 



10 



LU 




1 234 5 6789 10 

Vp 

FIG. 2: Bayesian evidence E(M) (solid red line) along with 
the unnormalized Bayesian evidence (green dotted line) the 
constant q model as a function of the prior volume V p . Nor- 
malization is arbitrary. 



E. Prior selection 

In general, prior selection is a crucial aspect of fitting 
model parameters and maybe even more when calculat- 
ing the Bayesian Evidence. From the definition, Eq. JJjJ, 
we see that the larger the volume of the prior parame- 
ter space, the smaller the Bayesian Evidence will be due 



The same quantitative behaviour is also found when 
considering more complicated models. It is hence clear 
that the prior volume and whether the maximum likeli- 
hood point is included within it, are critical when deter- 
mining the Bayesian Evidence of a model. If the prior 
is chosen in such a way the maximum likelihood point is 
not included, the evidence for that model may be orders 
of magnitude smaller than that of the same model with 
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somewhat larger priors. This may lead to incorrect con- 
clusions in comparing models, if one views the parameter 
priors as separate from the actual functional form of the 
model. 

In order to demonstrate how the conclusions may de- 
pend on the prior, we have plotted lni?oi for different 
prior volumes for the linear model in Fig. [3] From the 
plot we see that conclusions can slightly change with dif- 
ferent priors but radically different results can only arise 
if the prior of the linear model can be restricted to a very 
small volume. 
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FIG. 3: Bayes factor Boi between a constant and linear q 
models for different prior volumes of the linear model. The 
horizontal lines mark the significance levels according to the 
Jeffrey's scale. 



IV. CONCLUSIONS 



We have subjected the state-of-the art SNIa data to 
a Bayesian analysis, assuming an isotropic and homo- 
geneous universe, but making no assumptions about the 
theory of gravity or the matter-energy content of the uni- 
verse. The main results are that the two samples do not 
enable us to draw strong conclusions about the underly- 
ing model, but that there is no evidence that anything 
beyond a constant, negative deceleration parameter is 
required in order to describe the data. This is consis- 
tent with the conclusions drawn in Q based on a Princi- 
ple Components Analysis of the Gold sample. We think 
these results are important to bear in mind when assess- 
ing the relevance of fashions like 'phantom energy'. 
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